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ABSTRACT 


The thesis considers the practical prohlem of power system 
stabilizer design. 

An algorithm, which has not yet been used for design of 
power system stabilizers in the literature, has been used in 
this work. The algorithm involves partial and complete pole 
placement by output feedback. 

The algorithm leads to the design of a full rank static 
or a full rank, minimal order dynamic compensator and shifts 
the poles arbitrarily closed to a specified set of pole 
locations. It has been applied to a practical power system 
and static and dynamic compensators have been designed. 

The algorithm has been compared with the existing 
techniques for power system stabilizer design and advantages 
and disadvantages have been discussed, 

Itirthermore, the results of the two cases of partial 
and complete pole placement have been compared and merits and 
demerits have been discussed. 
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CHAPTER I 


IHTROEUCl'IOF 

The problem of improvement in dynamic stability is a 

fundamental problem in Control Engineering. In modem power 
systems, dynamic stability problem has arisen due to the fast 
acting excitation control systems. Though a number of attempts 
have been made in the design of -auxiliary controllers for 
excitation systems, known as 'Power System Stabilizers' and 
in many cases encouraging results have been obtained, • they 
are either based on heuristic ideas or on trial and error 
methods. In this thesis, a systematic procedure has been 
applied for the design of power system stabilizers. This ; 

procedure is based on a recent work in pole-placement using [ 

f 

[ 

output feedback [l0]. ' ' 

f 

The method used [jojis applicable to linear, time- \ 

invariant systems and is ’’full rank static or dynamic compen- [ 
sator design by output feedback” for achieving a given pole- 
placement. 

The thesis is organized as follows ; 

In Chapter 2 only those parts of Munro and Hirbod's work, ; 
which are relevant from the application point of view, have 
been given. The details and. proofs involved are given in part D 
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of Appendix A, Chapter 2 also presents a "brief review of 
existing pole-placement techniques. Chapter 5 introduces power 
system stabilizers and reviews briefly the existing methods for 
their design. In Chapter 4 is presented a number of power 
system stabilizer designs for a particular power system. The 
results obtained have been discussed and compared with existing 
results. Chapter 5 presents conclusions. 
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OHAPTER 2 

POlE-PEAOEiyiMT TECHFIQUES 

2.1 THE BASIC PRIRCIPIE 

The basic principle of pole-placement techniques is to 
shift the eigenvalues of the system to specified locations by 
state or output feedback. 

2.2 REVIEW OE DIE TECHNIQUES ^ 

In the case of state feedback it has been shovna that if the i 

[ 

system is state controllable, then all the eigenvalues can be 
shifted to arbitrary locations provided complex eigenvalues f 

exist in con jugate -pairs. If the system is not completely state ^ 

j 

controllable, then only those eigenvalues which are controllable, : 
can be shifted to arbitrary locations by state feedback pro- 
vided complex eigenvalues exist in conjugate pairs [2]. Many 
methods are available for pole-placement by state -feedback. 

However, in practice all the state variables may not be 
accessible. State reconstruction by state observers may be a 
solution of this problem. If '1’ represents the number of 
outputs, then the minimum order of the observer would be i 

(n-l), where n fepresents the number of state variables. * 

Outpht ; feedback without attempting state reconstruction is a 
more pradtfeai approach. ■ :E^en ' if a dynamic feedback is required, S 
its order will,! in' general, be less than (n-l). j 
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2.3 RBCENT TECHNIQUES 

Let us consider the linear, time-invariant system 

X = Ax + Bu 
y = Cx 

■where x is an n state vector, u is an m input vector and y is 
an 1 output vector, A,B,C are nxn, nxm and Ixn constant. „ . 
state, input and output matrices respectively. 

Lavison [1 ] has shown that if the system is completely 
state controllable and observable and matrix A is cyclic, then 
max(m,l) poles can be shifted arbitrarily close to a given set 
of specified locations under the constraint that complex poles ■ 

exist in conjugate pairs. The cyclicity of a matrix means that j 

its Jordon normal cannonical form should not contain more than 
one jordon block for any eigenvalue. But the cyclicity of A is 

caJiiL f 

not a strong constraint, firstly because of being a »a^^and | 

secondly because almost any arbitrary gain output feedback 
will make the system eigenvalues distinct and hence cyclic. 

This result is dp-e to Davison [3]. The basic idea behind the 
placement of eigenvalues is dyadic feedback .approach | 

which is discusaed .$n the next subsection. * 

2,3.d ' | 

■ ■ Ah' outlet' fhedba'Ck^^^^ is called dyadic feedback | 

if E is of the form ; . 



■where and are vectors of order m and 1 respectively, 

*tj 

Here, f^ represents transpose of f^. It is this dyadic 
structure of output feedback that leads to simple linear 
equatiom for solution if either, f or f is assumed apriori, 

CO 

transforming the multivariable system into either a single 
input, multi-output system or a single-output, multi-input 
system. Any direct attempt to evaluate mxl parameters 

of a non-dyadic F or (m+1-1 ) parameters of 

a dyadic P leads to the solution of non-linear equations -which, 
are quite involved even for simple cases. 

So the basic procedure is to choose f , if m < 1 or to 

w 

choose f if 1 m without affecting observability and contro- 
0 

liability and then to solve the n equations in the unknown 

elements of f or f ,as the case may be. Then the resulting 
0 0 

feedback matrix places max(m,l) poles to the pre-assign ed loca- 
tions. The procedure is given in detail in part B of Appendix A. 

The most obvious drawback of this static compensator is 

that rest of the poles ^ which are [ n-max(m,l)] in number get 

placed at random. So dynamic compensator comes into picture 
as a result of attempts made to overcome this problem. 

^ shown that for a single input, 
multi-output , syst^ compensator of order r * is 
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sufficient to place all (n+r) poles of the composite system, 
where *1^' is the ohservahility index of the system. They have 
also shown that a compensator of order r=(y'”1') is sufficient to 
place all (n+r) poles of the composite system in the case, of a 
multi-input, single output system, where 'V ' is controllability 
index of the system. 

Bras'Ch'and Pearson [5] have shown that a dynamic compensator 

of order r = min(VQ~‘', is sufficient to place all (n+r) 

poles of composite system to specified locations where 'V ' and 

'V * are controllability and observability indices respectively. 

The feedback used by Brascii and Pearson [5 ] is also of dyadic 

ihU<L 

type, but proof and method are different than of Ohen and 

Hsu [4] and more involved too. 

Subsequently, the work of Ahmeri and Vacro'ux [6] on dynamic 
compensators with dyadic output feedback shoi ’ j ed that 
max(^+r, ^+r) poles can be placed in an arbitrary set of pre"- 
assigned locations,’ where r is the order of the compens^ator 
and 

a =. Rank [B.AB, aS A^B] 

and p = Bank (A^)^c'‘, (a'*)^ 0*] 

Then* come the works using feedback as a sum of dyads 
and it. becomes ;'yideht that earlier works are only sufficient 
condition's fbr miiLii-vari&ble systems. Subsequent works oy 



Davison and Wang [7], Belletrutti [8], Munro and HirDod [10 ] 
show that more poles than max(m,l) can he assigned with static 
compensator and a lesser order of dynamic compensator than 
min(V -1, V^-1 ) place all of the (n+r) poles of the 

composite system, 

2,3,2 Two Dyad Approach 

Davison and Wang [7] have shown that with a sequence 

of two dyads min(n, m+1-1 ) poles can he placed arbitrarily 

close to a given set of pre-specified locations in ^almost all 

cases under the constraints that complex poles exist in 

conjugate pairs and system is completely state controllahle 
« 

,and ohservahle. 

The basic idea behind feedback by a sequence of two dyads 
is pole-retention as well as pole -placement. Through the first 
dyad, max(m,l) poles are placed just as in the sin^e' dyad 
case. Then during second dyadic feedback f or f is chosen 
to make [max(m,l)-1] poles either uncontrollable or unobser- 
vable, Then equations are solved for f^ or f of the second 
dyad to place as many more poles as possible. Davison and 
Wang [7] have shown that duri a g in ’almost all’ cases 
min(n,m+l-1 ) - [ max:(m,l)-1] poles can be placed during this 
second dyadic feedback. This leads to the placement ^ of 
min(n,m+l-1) poles in almost all, cases through a set of two 
dyads. But this static compensator still possess the old 
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dra-wback of placing the remaining n-(m+l-1) poles randomly 
if n > m+1-1 . 

To overcome this dra-whack dynamic compensators again come 
into picture hut vjith the nev approach of two dyad feedback. 

It turns out that a lower order compensator than the previous 
case of single dyad feedback is sufficient to place all the 
(n+r) poles of the composite system and this oMer r will lie 
in between 

[n-Cmn-l-D] , ^ 

where [,] represents the least integer, but greater than 
or equal to the enclosed quantity. 

Similarly in the case of partial pole placement of (q+r) 
poles [if n^q > m+1-1 ] of the composite system, where 
’r’ is the compensator order, it turns out that in 'almost all' 

cases y r will lie in between 

1^] ana [,-(..1-1)] , 

where [ . ] represents least integer, but greater than or 
equal to the enclosed quantity. But in this case of partial 
pole placement the- unassigned (n-q) pel -vs will, get 
placed randomly. 
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2.3.3 Extension of Two Dyad Approach 

Eurther extension of the two dyad sequence approach has 
heen made [8-10] hut increase in the nuniher of dyads beyond two 
has not resulted in lesser order dynamic compensator for 
complete pole placement. Similarly, increase in the number of 
dyads beyond two has not resulted in assigament of more poles 
than min(n,m+l~1) for the static compensator case. 

But Belletrutti [8] has recently shown that increase in 
the rank of feedback compensator results in lowering of eigen- 
value sensitivity with respect to parameter variations. This 
has been confirmed by Gomathi [9] through a numerical example, 

2.3.4 Eull Rank Compensator by a Minimal Sequence of Dyads 


Munro and Hirbod [10 ] have proposed a new algorithm for 
full rank compensator design by a minimal sequence of dyadic 
feedbacks. In essence, it is an extension of pole-retention and 
pole-placement technique. It possess all the advantages of 
earlier techniques besides its own. Since every feedback is 
of dyadic form, it preserves 'the mathematical simplicity of 
linear equations. Since it is of full rank, eigenvalue 
sensitivity, due to parameter variations, is less. 


: <xr\ axkmfjiyn 

furthermore, since the feedlxclifi^; is comp o s ed .of ma x im al 


ef dyads, 


ae-qumee &£ dyads, so it result^ minimal order dynamic 

’,"V ' '' ■ 

I 

compensator^ 8^ -fe aii d^-no-’yy^ 

IwifvXck avvA Jjd'cjvyyU'nJud a, d^&jdjUs . 
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The hasic idea in Munro and Hirbod’ s algorithm for the 
case of m 1 is given helow. 

In the first step is chosen as [l ,0,0, ...» and 

(1 ) ° 

f^ is found to place as many poles as possible. 

In the second step f^^^^ = ff ^ ^,f 2 ^, 0 , . . . jOj"^ ; 0 is 

chosen io make maximum of the assigned poles in first step 

(2) 

uncontrollable. Then f^ is found to place as many more 
poles as possible to specified locations. 


(i) 


This process is repeated until the (m-1 )th time. Each time 
“ •••» * »0 , , 0 J^ji •••» (i3i~1)jf^' ^ 0; 


is chosen to make maximum of the previously assigned poles 
uncontrollable and then is found to place as many more 

poles as possible to specified locations. 

(m) m „ m 


In the mth step •••» ^ 

is chosen to retain (m-1) previously assigned poles and then 
f is found to place the remaining poles. 


m 


If number of poles to be assigned/5 1-4: (m+1-1), then a 
constant -will result and E -will be static. If n>q>m+l-1, 

the resulting will be dynamic and E will also become 

dynamic. 


This method is given ih detail in pait D of Appendix A, 
Eor the case of m T' 1, its du^l is applicable. 
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This method is quite suitable lor partial as well as 
complete pole-placement and in each case the method results 
in a minimal order, full rank compensator. In case of 

^ m+1-1, the resulting compensator is static in' almost all' 
cases. 

Though in the case of partial pole placement the unspeci- 
fied poles get placed randomly, it may so happen that they may get 
acceptable locations. In such cases, partial pole-placement 
is at an advantage in terms of lower order compensator but with 
good system performance. 

But the cases, in which locations of unspecified poles 
are not acceptable, complete pole-placement can always be done 
for proper pole-placement. However, possibility of obtaining a 
good design with partial pole-placement should always be explored 
first , 

In the present thesis, design of compensators have been 
done for both the cases, namely, partial and complete pole 
placement. In many cases of partial pole-placement results are 
very good and in a few cases even static compensators have 
resulted. But in a few cases poles have got placed at odd 
locations too. The complete pole-placement has always resulted 
in satisfactory results. ' 
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CHii'TER 5 

PO¥M SISTH4 STiBILIZERS 


3.1 INTRODUCTION 

Attempts have been made to improve the transient 
stability of power systems by incorporating fast acting, 
high gain excitation systems with synchronous machines. 

Althoi;}gh this leads to thre improvement in transient stahility 
of the system, it often gives rise to the problem of oscilla- 
tions and dynamic instability depending upon loading condi- 
tions [17-19]. Oscillations of sustained or growing nature 
have been reported in the literature [18-21]. The basic 
cause of such oscillations is the high gain of the voltage 
regulator, which alongwith the fast acting static- exciter 
introduces negative damping in the system affecting its 
gyna'^ic stability [l6]. So a need for another supple- 

mentary unit, to be used in the excitation control arises 
which can improve the dynamic stability of the system [17]. 

This supplementary unit is termed as Power System Stabilizer 
(PSS). 

PSS is an auxiliary controller which receives a feedback 
signal from rotor angular position, angular velocity, frequency 
or acceleration and provides corrective signals at the input 
of the excitation system in order to damp out the oscillations 
in the syst^^ 
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Since angular velocity, rotor angle, acceleration and 
frequency are directly measurable variables, only these are 
used as feedback signals. 

Since dynamic stability of a power system is related to 
the changes in the system response due to small disturbance in 
the system, the PSS can be designed with respect to the 
linearized model of power system about an operating point. 
Clearly, pole-placement techniques using output feedback can 
be used for PSS design. 

3.2 PSS DESIGN OBJECTIVE 

As discussed in Section 3.1, the reason of instability in 
power system is negative damping introduced by the high gain 
excitation control system for system voltage regulation. 

Since in a linear system negative damping can be viewed 
as being due to positive real part of a complex eigenvalue, the 
objective of PSS design will be to shift the unstable complex 
eigenvalues to appropriate locations in left half plane. This 
will make the system stable. 

3.3 A BRIEF SURVET OP EXISTING PSS DESIGN TECHNIQUES 

Initially classical control theory [l5,l6] was applied to 
get improved system response. But it was limited to single- 
input, single- out put systems. 

Later on, modern control theory was applied. Here two 
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main approaciies are available, namely, that based on lihear 
optimal regulator theory [ 22 , 23 ] and that based on- pole- 
placement techniques [9,24,25]- 

Application of linear optimal regulator theory results 
in a sta.te feedback. However, a drawback of this approach is 
that all states have to be accessible for measurement. Further- 
more it is an iterative technique in which one has to keep on 
changing the weighting matrices Q and R in the performance 
indez tiJ-l satisfactory response is obtained. 

Among the pole-placement techniques, state feedback has 
been used for PSS design [24,25,28]. But it is not a practical 
case, because all of the state variables are not available for 
measurement in power systems. 

Among output feedback pole-placement techniques [9], 
there are two well known approaches, namely, algebraic output 
feedback and dynamic output feedback. 

In the case of algebraic feedback, a dynamic compensator 
with defined structure and parameter is chosen first. The 
choice of compensator structure and parameters usually depeild 
on experience and heuristic ideas. After choosing this compen- 
sator the algebrliic feedback law is found to place 

the poles to appropriate locations [9]. 

In the case of dynamic feedback, too, a compensator 
structure of some fixed order is assumed first and then the 
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parameters of this compensator are calculated, to place the 
poles to appropriate locations ['9]. 

In all of the above methods some drawbacks are always 

OrfX 

present. They are either based on heuristic ideas or^trial 
and error methods or they suffer from the preassumption on 
compensator structure and order or from inaccessibilily of 
states, ho straightforward method, which can directly give 
the order, the structure and the parameter values of PSS has 
been used so far. 

' In this thesis su.ch a method has been applied for PSS 

design. The algorithm used in the method is due to Munro and 

Hirbod [l0]. The algorithm is applicable to partial as well 

iy\ rfc 

as complete pole-^p la cement a par % ^ fro m thre advantages men- 
tioned above. In the case of partial pole-placement the 
compensator will turn out a static one if the number of poles 
to be assigned is less than or equal to (m+l-l), where m and 
1 are the n\imber of inputs and outputs respectively. 
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CHAE'TER 4 

POWER SYSTEM STABILIZER DESIGN 

4.1 INTRODUOTION 

It has already been shown in Chapter 3 that PSS can be designs 
with respect to the linearized model of power system, at the 
operating point of interest, for improvement in dynamic 
stability. So in Section 4.2, a linearized state-space model 
of a synchronous generator connected to an infinite bus through 
a transmission line has been chosen. Derivation of this 
linearized model is given in Appendix B. In Section 4.3, 

i, 

the problem of PSS design has been formulated in precise f 

terms. In Section 4.4, the algorithm to be used for PSS j 

design is given. In Section 4.5 numerical results have been | 

t 

presented and in Section 4.6 the results have been discussed. ■ 

I 

4.2 POWER SYSTEM MODEL [ 

The power system, whose dynamic behavioixr is to be 
improved by output f eedbach,. copsists of a synchronous, .genera- 
tor connected to infinite bus thtpugh a transmiesion line alongwi*; 
excitexrvolta^e' r^gnintor as given in ^ Pigs* 4.1 4»3, , 

Appendix B. gives ,, the , nonlinear equations governing the system ^ 

dynamics ;,;'ani/;tHe,.l|i:^ari^ei ' equations de^^^ the non- . 

lin ear equations ' ah Out an Op erating p o int . The lin eari z ed i 

sy et em . equatio.ne; .are-,. ' 



f $S ■ 

f. 
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X = ■+• Bv 

y = Cx 

"tr 

where x = ^w, a,x^ ] and v = [ 4. u] . Here x is an 

n state vector and v is a scalar input, y is a scalar/vector 
output depending upon the number ot outputs. 1,3,0 are 
constant matrices of appropriate dimensions. A and B are 
given below % 

I "■ 0.16312 0.3773 0.0 - 0.34722 

I 0.0 0.0 1.0 0.0 

A = ! 

- 9.0999 - 25.7124 0.0 0.0 

2983.2839 - 2436.8389 0.0 - 100.0 


0.34722 

0.0 



There are three different cases of output feedback which 
have been examined. They differ in the choice of output 
variables. The three cases are as follows j 

a) feedback from angular velocity (hw) only 

b) feedback from angular velocity (aw) and angular 

position ( A 5), and 

c) feedback from angular acceleration ( Aw) and angular 

velocity ( A w) . 
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Only the ahovc-ment ioned variables are considered for feedback, 
since only these can be directly and easily measured. In all 
the three cases matrix C and output y are different and are 
given below ; 

Case (a) C = [0,0, 0.0, 1.0, 0.0] and y = [y^] = [ 4 w] 

[O.O 0.0 1.0 o.o“ 

Cas e ( b ) 0 = and y = 

0.0 1.0 0.0 0.0 




placement have been attempted. Eigenvalues of the system 
without PSS is given in Table 4.1. It is seen that dominant 


eigenvalues ai-e unstable while the remaining two are stable. ? 

r 





So the main, design objective while designing PSS is to 
shift the real part of the dominant eigenvalues to a location 
in negative ha2.f plane such that system with PSS becomes 
stab3.e. 

The new locations for the dominant poles have been chosen 
as -1.0 + j5.5878. So, the minimum number of poles to be 
placed is 2. Now different possible cases of partial and 
complete pole-placement are given below ; 

Case I ; Partial pole placement ; Since n = 4, hence 
different possible cases are (a) q = 2, and 
(b) q = 3, where q poles are to be placed. 

Case ^ specified locations are given below. 

(i) -1 .0 + j5.5878 

(ii) -1.0 - j5.5878 

Case (b) q = 3. The specified pole locations are given 

below. 

(i) -1 .0 + 35.5878 

(ii) -1 .0 - 35.5878 
(lii) -12.46339 

Case II ; Complete pole-placement s In this case q = n = 4, 
and the specified pole locations are given below ; 

(ii) V-tvO 3:5.5878 
: (iii): 

, (iv) -88. 226,2 
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4 . 3 PROBLEM E ORMUIATI ON 

There are three cases and in each case the prohlem 
formulation is different. These cases arise due to the 
partial or complete pole-placement and due to static or 
dynamic type of compensator, let q be the number of poles 
to be placed and m and 1 be the number of inputs and outputs 
respectively. Then in the present power system case 
n > m+1-1 . 

Oase I s q m+1-1 ; partial pole-placement with static 
comp ensator 

The problem is to find an output feedback law 

u =-E,y 

where P is a constant matris , such that q poles of the 
original system get placed at the specified locations. 

Oase II ; n > q, > m+1-1 s partial pole-placement with dynamic 

compensator 

The problem is to find an output feedback law 
U(S) =:-P(S) Y(S) 

where F(S) is a proper rational polynomial matrix of minimal 
order r, such that (q+r) poles of the composite system get 
placed at the specified locations. 
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Case III ; n = iii+l~1 •; complete pol emplacement with. 
dynamic comTDensator 

The problem is to find an output feedback law 
U(s) =-T'(s)-T(s) 

where P(s) is a proper rational polynomial matrix of minimal 
order r, such that all the (n+r) poles of the composite system 
get placed at specified locations. 

4.4 ALGORITHM 

To determine a P or P(s) as required in Section 4.5 > the 
algorithm of Munro and Hirbod has been used. Necessary details 
and proofs are given in part D of Appendix A. 

4.4.1 Preliminary 

This algorithm leads to the determination of a non-unique 
but full rank feedback matrix P and places the eigenvalues of 
the closed-loop system arbitrarily close to the specified 
locations. Pirst of all the of n^sm+l-l has been considered, 

The desired matrix P can be constructed by a minimal 
sequence of dyads, that is, 

^ (i) f U) 


¥= A 




0 


*• 4.ct An 


where t A* v eo 


!ir#qt6rVfu^^^ is an 1x1 vector 


and 


H s mih(m,l) The ipatrlx P can equally be expressed as a 
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product of two matrices 


F = P F 
0 0 

where columns of F are the vectors f and row of F^ are 

( ■ ) ° ° ° 
vectors f^ »for i = (1,2, p) in "both cases. ITow the 

algorithm is given helow for the case m 4 1. For 1< m, the 

dual is applicable. 


4*4.2 Algorithm for the case m4 1 


Following are the m stages in the algorithm. 

(1) Select f as [l,0, ..., 0]"^, and determine the 

corresponding to place as many pole as possible 

in desired locations. 


(2) For £2^' ***» determine f^^ and 

f„^; f^^ ^ 0 to retain as many previously assigned poles 

(2) 

as possible and determine f^^ ‘ to place at least one 
more pole in a desired location. 

(3) For = [fi^» £ 2 ^’ £^^,0,..., 0]^, determine 

i "^2 ^ £3^; £5^ 0, to retain as many previously 

( 3 ) 

assigned poles as possible and determine to place 

at least one more pole in a desired location. 

Repeat the sequence till f ^ and are determined. 


(yn) 


For 

■ ■ o . 



f ®], determine £.^ to £ 
la 1 


f 5^ 0, to retain (m-*1 ) previously assigned poles and then 

m ■ f ■ ■■ 
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determine to place min(l,n~m+l) poles in desired 

locations which, is possible in 'almost all' cases, [7] • 

So all the n poles have been placed for the case 

q n ^ m+l~1 . In case of n q > m+1-1 , the first (m~1 ) 

stages are carried out in accordance with the sequence given 

above and in the mth stage f is determined as usual and 

c 

then all the remaining poles (yet not placed) are placed in 
according with the method of Chen and Hsu [4], which is given 
in part 1 of Appendix A. 

4.5 B.ESULTS 

The results are given in Tables 4.3 " 4. 5 • 

Unit step responses of rotor angle and angular velocity 
with respect to time have been given in Pigs. 4.4 and 4.5 
respectively for all causes except two. These two cases 
correspond to those partial pole-placement cases in which 
the unassigned poles have got placed in positive half plane 
leading to instability. The various controller configurations 
for which responses have been plotted are given in Table 4.2. 

Different parameters corresponding to various cases have 
been defined below : 

Case (a) s The .cas e, bf partial pole-placement ; 

Here input-out “put relati^o^^ 
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Table 4*2 

Various Controller Configurations 


Type of compensator and feedback variables 



1 . Static compensator; feedback from cw and 

2. Static compensator; feedback from aw and 

3. Dynamic compensator of order 1/ Feedback from aw only 

4. Dynamic compensator of order^^l; feedback from 

AW and aw 

5. Dynamic compensator of order 2; Feedback from 

AW and Aw 

6. Dynamic compensator of order 2; feedback from 

A w and A ■'S 

7. Dynamic compensator of order 3; feedback from aw only 


where f == and f^ and fg are constants. 

Case (b) : The case of complete pole-placement : 

(l) The case of 1 = 1 j In this case input-output 
relationship is as follows 

tr(s) =-f(s) ^ T(s) 

1^0 + Pi 

where f(s; .= ' 

•„ S ' 't ^'1 "t 

where r is thb *5rd<^? b'f cbmpois^'^^t’* . for the first order 
comp ehsater p ; ito^iK' 
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^ s 4- [3 . 

F(s) = -2 ."1 , 


s + oc 


for the second order compensator form, it is of the form 

^2 

P( s) = 


t^oS + Fi^s + ^2 


s + a ^ s 4* q:^ 

and for the third order compensator, it is of the form 

'(i + p s + 

F(s) 1 ..2_^ 

s + a. s + a^s + 

1 2 5 

(ii) The case of 1 = 2 s The input~output relationship 


U(s) =-F(s)-T(s) 


where F(s) is of the form 


^'(s) = [ r^(s). ipCs)] 


with ^ ( s ) 


and 


r r-'l r“*2 

s + a^s + a^s + . . . + 

Pol®"' + Pll®"’”' + + ••• + Pri 


for i = 1 f 2. 


There are two cases of r and they are as follows, 
For r s=, 1 


F(s) 


1 


[Pm® + Pii» Po2® ^ ^2^ 


¥+aj ‘■'"or ■ '"11’ ro2‘ 
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Ta*le 4.3 

btatic Compensator (Partial Pole Placement ) 


S. fv^. 

1 


Feedback 

variables 

AW & 

/^^w & Aw 

No. of poles 
to be assigned 

(q) 

2 

2 

Assigned 

locations 

-1 .0+05.5878 

-1 .0+05.5878 

No. of poles 
left for random 
placement (n-q) 

2 

2 

Valu es of the 

compensator 

parameters 

f^= -7.60 

f^= -8.83 


f2= 19.84 

t^= -0.62 

Location of the 
assigned 
poles in closed 
loop system 

-1 .000 +05.5878 

-1.000+05.5878 

Location of the 
unassigned poles 
in closed loop 
system 

-88.28, -9.894 

-87.97,-12.14 

_ _ _j 



S./Vo. 

o. of poles 
0 lie placed 
q+r) 

re -as si gned 
ocati ons 


eedlack 

arialles 

ompensator 

arameters 


ocation of 
,s signed 
oles 


location of 
nas signed 
oles 


Talle 4 *4 

Dynamic Compensator (Partial Pole Placement) 


... __ 

r^i 

^+r=3 

^•C>135.5878, 

31.46339 


A ¥ 

^ .,=-855.01 
o( =99.092406 


-1 .000001 
±35.5877, 
-12.46339 


-85.55061 , 
-99.24153 


c3> 

q=3 
r=1 
q+r =4 

-1 .0+35.5878 
-12.46339, 
- 88.2262 


A'w &, aS 

^01=°-°' 

=1 1 . 24 

2=278.06 
o(^ =1.5669 

-0.9967508 
+ 35.60701 , 
-1 2.46104, 
-88.22741 


0.9518578 


•3 

q=3 

r=1 

q+r=4' 

-1 .0+35.5878 
-1 2.46339, 
- 88.2262 


A w& Aw 

/^o2~® 

^=-37.27 

/ 1 ^ 2 =- 407.05 
^ 7' 01^ 

-1 .000001 
+35.5878, 

- 12 . 46339 , 
- 88.22 6 2 


-45.644459 


h 

q=3 

r=2 

q+r=5 

-1 .0 + 35.5878 
-12.46339, 
- 88.2262 
- 20.00 

Aw 

/^o=6.38 

/3^ =624. 24 
^ 2 = 7013.42 
^;,= -15.2316 ; 

^P= -777.6149 i 

-0.999999, I 
+35.58779, i 
-12.46363, 
-88.22621 , 
-19.99979 


37.75814 


Nofa !- C6m^e.7)Aado^ o^ckit. 
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Table 4.5 

Dynamic C ompe ns at Complete Pole Placement) 



i 

^ 

; .._3. 

"'Feedback : 

variables . ; 

Aw &aS i 

Aw &A<^ ; 

^Aw only 

No,,, of poles i 

n=q=4 i 

n= =4 1 

n=q=4 

to., be i placed > 


• r =2 [ 

- r=3 

q=n,r,. n+r ; 

n+r =6 1 

n+r =6 ; 

n+r=7 

Pre assigned 1 

-1.0+35.5878 I 

-1.0+j5.5878 , 

-1 .0 +.3 5''. 5 87 8 , 

locatlbns ; 

-12.46359 1 

-1 2 . ..4.6339, : 

-1 2 . 46339 , 

-88. 2-2 62c 1 

-88.22 62.': i 

-88^2262/ 

t 

i 

-2,0.00 fe-l 5.00' 

—20.00 .& -1 5 .00 i 

“»**20 *^00 


' * ' • 

• ^ ^ i ’ ♦ i 

--i- 9 1.00 & - 25.00 

Order of 

r =2 i 

r =2 ' : 

r=3 

compen- ' 




sat or. ' 




Compensator , 

^oi=-2.46, i 

Boi=-0.23, : 

M=0.0, ■ 

p= - 8.32 

parameters ; 

^^ 2 =^. 0 ,’ 

A -=-240.67, 

p -88.60 


C. 

,*11=2,5.92, 

% 2 = 0^0 * ’ 

p^^=-493.07, 

:^^=-^572:79. 

■^ = - 14297.32 


r \ \ . . r C 

,^i2=o.o;'- 

V^1=-19'70.51, 

B22=®^57'.78' ’ 

^_= L 5 6262 164 
'= 62.52647 
% = 1302.5721 

. ^ - 1 ’■'+ i 


=.-37.52647, 

'=36.786575, 

"ci ,=308 ; 91 7 

\ 4il 

pC^ = 7722.9402 


=38^:967'34' 

^ ’’’'ll 4 ' f i ' 


Locations of 
aesigaed 
poles,' .0 

-i .000002 

+t)505878;! .. 
-1;2.4634tv 
188.22621, 
-20.00044 , 

-T-0±d5.5878, 

-■■-12wi4633-9V7o. 
-88 . 2262/ : 
-49.99999, & 
-15.0006t> 

-0.9999847 
-+35'.587785, 
-r2.4651'1 , ' 
-.= 88.226 2/.^ 

— 1=9 9 9 88 6 , 

A c: fic\ A rv4 > ^ 



. _,-1 t; , {J.'SCIO 1 

3 aUU 1 U r ^ 


-1 4»yyy^ ^ 

...PAtVcO'P,' '; ' 


«2f # 00045 ^ ■ 


'1 " - , / , ,, , 


4-,....-. 







troUer conf. t, 3. 4,5^ 7. 
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Step respond |ar aT)gulqr velocHy,£iW, 
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and for r = 2 


P(s) 


3 H-a^s+a^ 


^^ 01 ® ■*■ ^ 21 ^ ^^ 02 ® ^ 22 ^ 


The values of different parameters of P as well as P(s) are 
given alongwith in the tables 4.2 to 4.5. 


4.5 DISCUSSION 


After a careful inspection of the results_s it appears that 
w is the most effective feedback control element among all 
the three measurable variables, namely, A6, £-w and aw. 

Purthermore, it also appears that a combination of aw and 
n w as feedback elements is a best choice. 


The case of partial pole-placement reveals that the 
unspecified poles can get placed anywhere which is in accordance 
with our expectation. But in most of the cases they have got 
placed at good locations and in a few cases at extremely good 
locations, as in the case of Aw feedback for placement of 2 
poles, and in the case of aw and aw feedback for placement of 
3 poles. 

Static compensators have resulted in the cases of a w 
and 4 6 feedback as well as Aw and aw feedback for placement 
of two poles and the unspecified poles have also got placed 
at reasonably good locations. 

Hence attempts for partial pole-placement are not fruitless. 
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As far as tiie complete pole-placement is concerned, it is 
found that the pole-placement is more or less exact, further- 
more in tho complete pole-placement, case there is a certainiy 
about the improvement in system dynamic response because there 
is no chance of getting any pple plac^ed at an tindesirable loca- 
tion. 

■"As far as the algorithm' is donberned, it is simple and 
straightforward and places, the '.poles quite accurately, further 
more, this algorithm itself determines the minimal order 
and structure of the compensator which is the most important 
advantage of this algorithm over the other existing ones. 

A study of the responses given in figs. 4.4 and 4.5 reveals 
that there is no distinct advantage of the dynamic compensators 
over the static ones, because the responses in all the cases 
are almost identical. In the present case it is so because 
new poles and zeros, originated due to dynamic type of com- 
pensator, are too close and cancel the effect of each-other. 

But this type of behaviour is not guaranteed in all cases. In 
general, all the three cases, mentioned below can take place. 

(a) JiTew zero is nearer to doaiinant pole than the new pole. 

(b) lew zero and new pole are nearly equidistance from the 
dominant pole. 

(c) Hew zero is at greater distance with respect to new pole 
from the dominant pole. 
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Since^ the dominant poles are placed at the same loca- 
tions in dynamic as well as static compensators and original 
zeros do not shift hy output feedback in the case of single- 
input, multi-output or multi-input, single-output system. The 
response of a dynamic compensator can differ 
that of static compensator only due to the effect of these 
newly added pole-zero locations. 

Suppose new zero and pole are at z^ and p^ respectively 
and the dominant pole is at p^. Then the response corresponding 
to the most dominant pole of the system will be 


^d ^n 

Pd “ Pn 



where K e is the response 
pole of the^ static compensator. 


f 


due ' to' the dominant 
Hence the response of the 


dynamic compensator will be better than static one, if 


z 


n 


^ P 


< 1 . 


n 


So, only in case (a), the dynamic compensator will 
result in better response than static one. In case (o) the 
response will be nearly identical, which is the present case 
too and in case (c) response of dynamic compensator will be 

worwse than the static one. 

Hence, dynamic 'compensator will be advantageous over 
static one if new zero is nearer to dominant pole than the new 
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pole. But ;.^.v it oan not be guaranteed that the new zero 
will always be nearer to dominant pole than the new pole. 

Hence it can be concluded that 5 ii jiot 
necessaiy that the dynamic compensator will give better res- 
ponse than the static compensator if the unspecified poles 
of the latter got placed at acceptable locations. 

So, the superiority of the dynamic compen- 
sator with complete pole-placement over lower order dynamic 
or static compensator can be claimed in only those cases 
where the unspecified poles of latter system do not get 
placed at acceptable locations. 
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CHAPTER 5 
OOFCEUSIOlf 


In the present work a new algoritiim liss been used for 
PSS design. This algorithm is better than those which have 
been used earlier for PSS design in many respects, besides 
being simple, straightforward and computationally easy. It 
is neither based on any heuristic procedure ' nor involves 
trial and error procedures. Ho apriori assumptions need bo 
made about compensator structure. Apart from these, the 
algorithm directly gives the minimum possible order and 
structuz-e of compensator whereas earlier works , dealing with 
PSS design, involve assumptions regarding the order and 
structure of the compensator.. 

Besides, presenting a new algorithm to the PSS design 
area, some other conclusions have also been reached regarding 
the suitability of a dynamic compensator.- It has been con- 
cluded through .l ogic end example that it is quite fruitful 
to attempt partial pole-placement and it is quite likely that 
the unspecified pole may get placed at acceptable locations. 

It has also been shown that if a lower order dynamic 
compensator (due to partial pole-placement) exists, then it is 
not necessary that the system with higher order dynamic 
compensator (due to complete pole-placement) will always have 
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better response than that with lower order. In general the 
response may worsen, may remain same or may improve. It all 
depends upon the newly added poles and zeros. 

It has also been concluded that dynamic compensator xfith 
complete pole^-p la cement is better than the lower order 
dynamic or static compensators with incomplete pole-placement 
only if the unspecified poles of the latter cases do not get 

placed at acceptable* locations . 
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iPPBroiX A 

VARIOUS OUTPUT FEEDBACK APPROACHES 


GEHERAE 

Let a linear, time-invariant system be given by 
x = .Ax + Bu 
y = Cx 

where x is an n~state vector, u is an m~ input vector and y 

I 

is an 1-output vector and A,B,C are constant matrices of 
order nxn, nxm and Ixn respectively. 

How this system can be represented in the frequency 
domain by the transfer function G(s) as 


G(s) - C(sI-A)“'^B 


C ido(sI-A)B 
det(sI-A) 


r 


/ffe( 3 


S 


and 

Y(s) =:G(s)U(s). 


where = det (sI-A) is the characteristic polynomial of 
the system and its zeros *•*> '^n^ poles of 
the system. 


How our problem is to determine a feedback law of the form 
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a(s) = R(s) - R(s: Y(s) 

which assigns the (ii+r) poles of the composite system to 
specified locations 

Here, Il(s) is a vector of r reference inputs and r is the 
order cf I'(s), the compensator. 

How it is well-lcnown that any direct attempt to evaluate 

the parameters of I'(s) results in a set of non-linear equations 

It is equally well known that if I is dyadic, that is, P is 

"t? 

defined as the outer product of two vectors f and f^ as 

CO 

1 ? = 

in 1 

where f C then a direct attempt to solve (m+1— 1 ) 

c o 

independent parameters defining P also requires the solution 

of a set of non-linear equations. The usual approach to 

this latter problem, is to choose either f^ or without 

affecting controllability and observability. It is this 

"tl 

choice of f or f that makes the remaining part of the 
system design linear. 


PiPT A 

DYIDIO COMPEHSATOR DBSIGH ( SIHGIE-IHPUT , MULTI-OUTPUT 

CASE) [4] 

Consider a system described by an irreaucible 1x1 strictly 
proper rational transfer function matrix G-(e), defined as 



= A^) 


i N^(s) 


%(s) 


^oC^) 

Ar,(sT 


where A^(s) = s^+a,s^'‘^ 


+ - . . + a 


n 


and 


Bi(s) = Pi, s”-' + Pi2s“"2 + 


'*’ ^XT) » (i - 


1) 


is the least common denominator of; vthe elements of Gr(s), 
and , theref or-e, is alsothe characteristic^ polynomial of Gr(s). 

problem is to determine a ix:l dynamic feedback compen- 
sator F(s) having least degree such that the closed-loop system 


H(s) = [I 

has a desired set of 
Now by assuming 


+ &(s} F(s) G-(s) 

poles (?Vjr ; i = 1,2, 
F(s) of the form 


n+r) . 


I'(s) = 


1 


A. ( s ) 

c'' ^ 


[M,(s), 


, M3^(s)] = 




where 

A^(s) = s^+h^ + . . , + h^ 

and 

' »i(s) = «!<, + Si, + «12 + ••• + Sir 

jL ••#jr 1# 


the resulting closed loop polynomial A^(s), defined as 
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Aj(a) = ... + ... + d„^^ 

can be written as 

1 

A(i(s) = ^^(s) A^(s) + 2 %(s) 

i=1 

The above result can be proved as follows ; 

r^cs) 


H(s) 


[I + G(s) PCs)]"^ G(s) 


G(s) [I+FCs) G(s)]~^ 


(A1) 


_ G(s) M.in + P(s) G(s)' 

_ -.(iet [l+PU) G(s)J 

But F(s) G(s) is a 1x1 matrix, hence a simple rational poly- 
nomial and it is well known that adjoint of a single element 
matrix is 1 . So, 


H(s) 




G(s) G^s ) 


^ ( T = ditlH-Pfs) G(s)J == 1+F(s) G(s) 
("(s) 

"^TsT A(.(s) fo^s) 


1 + 


’Aq(s)Ao(s) 


Hence, 


A^(s) = A^(s) Avq^®) + FqCs) r^(s) =A^(s)Aq(s) + 

1 

^ M^(s) l^(s) 
i=1 
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Now by equating the coefficients of powers of s of eqn.(Al), 
a set of (n+r) equations in [(r+1 )(1+1 )~1 ] parameters of F(s) 
will be obtained xfhich is given below : 

Z P = 6 (i2) 

where the vector P contains the parameters h.'s and ^ ' s of 
P(s) and the difference vector 6^ contains the coefficients 
6 . of the polynomial 

J 


6^(s) = A^(s) - A^(s) s^. 

Now, X^, P^ and 6^ are given below : 
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Now, for eq[n.(A2) to have solution, 6^ must lie in column 
space of X^, i.e. 

rank(X^,6^) =:rank(X^). (13) 

But, this condition will he satisfied, if the number of 
Independent columns in X^ will be more or equal to (n+r), the 
number of rows. 

From here the minimum possible value of r can be calculated. 
The condition for this will be 

[(1+1 )(r+1 )-l] ^ n+r . 

On simplifying, 

> [Sji] (A4) 


r 
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where [K] implies the least integer ^ K. 

Similarly the upper limit can also le found with help of 
the results of Ghen and Hsu [4] and Brash and Pear's® on [5j and 
it turns out that 

r ^ [n~l] 

Hence r is bounded as 

< r [n-1], (A5) 

where [K] implies the least int^er E. 

Here, one point should be noted that condition (A3) may 
get satisfied for r<^ special cases Q'f 6^. 

Hence, the best way to find minimina value of r is to 
follow the method given below ; 

a) Set r = 0 

b) Check condition (A3). If satisfied, go to stejp d. 

c) r <^r+1 • G-o to step b. i 

"'min = I 

After determining r by above method a solution for equation I 

(A2) can be given by 

■ . f' 

and [ 

f’y "■ > ' i 

1 

E 

F 
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For the first case solution is unique, while for the second 
case it is one of the non~unique solutions. 

PART B 

DYADIC C0I4PMSAT0R DESIGN (MULTIVARIABLE CASB)[ l] 

Here the procedure is to choose f^"^ ; if 1 < m or f^ ; 
if m 1 (see equation Gl), such that the resulting system is 
completely ohservable and controllable. The corresponding 
matrix having dimension n x inaz(m,l)j,will have full column 
rank max(m,l) and 

I'd = ^0^"’ C 

will give the best solution for in least error square 
sense. Furthermore, if 6^ lies in the column space of X^, 
then the solution will be exact and unique. 

Now, let q n be the number of closed-” loop poles to be 
assigned, while t = n~q poles are allowed to assume arbitrary 
values. The closed-loop characteristic polynomial s) for 
the present case can be written as 

A^(s) = A^(s) A^(s) 

= (s^+d^s‘^~Vd2S^~^+. . .-i-dg^)(s'*^'+e^s"^"'V. . .+e^) 

(B1 ) 

and the difference vector dgin the present case will be 


redefined as 
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/- 


5o- 5^ + + e.p6^ f ... + 


where 


and 


6q(s) = ^^•^(s) ' s 


2 2 


A^(s) 




(B2) 


6^(s) = Aq^(s)- s 


t”-i 


Then the eq^aation (A2) becomes 


i = 1 , . . . , t , 


^ 0^0 = «o + D e 


(B3) 


where 


and 




f e^ > ®2 ^ » ^t ^ 


Here vector e"^ is transpose of vector e. Eqn, (B3) can fur- 
ther be manipulated to inr the form 


A '■V A 

^0 = «c 


where 


p.. 


=-- [S^,-D] and 




P 


0 ! 


(B4) 


Here dimension of is n» [mas:)m,l)+t] = n*[n~q+max(m,l) ] . 
How one of the following three cases may arise ; 


Case I ; If q = max(m,l), then the dimension of is 


n*(q+t) = nxn, and can be uniquely determined as 

A A 1 A 

Pq = [Xq] ^ 


(B5) 
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Case I I : If q < max(m,l), tJien number of columns 
[n+max(m, l)~q] will be more than nxunber of rows and a number 
of solutions are possible. One of them is 


A. _ A + A _ 1 A 4- ^ 

X^] Xq^ 6^ 


(B6) 


Case III % If q ^ max(m,l), then the number of columns 

[iir*{q+max(m, 1) I ] will be less than the number of rows, and a 

A 

solution will exist, if 6^ will lie in the column space of 

A 

X|^ and solution will be 


P 


0 


A ^ A 

= 


Vr’ 


(B7) 


A A 

If does not lie in the column space of X^, then a 
dynamic compensator will be needed. In such a case, 6^ of 

A 

eqn. (A2) should be redefined to give 6^ as 

6^(s) = s"*^ - s^ (B8) 


and in this case t' = n+r-q. Under such condition equation , ' 

(B4) will become 

[' 

A A A ■ ' 

X P = 6 (B9) 

r r r ; 

A ^ 

The dimension of matrix X^ will now be (n+r)- [max(m,l)* (r-i-1 )+r+ 1 [ 

Upper and lower 'bounds in the present case can be deter- | 

mined in the same way, as given in part A . It turns out that j 

c i' 

for placement of (q+r) poles the compensator order r is bounded ^ 


as 




A— 1 1 



r .<1 [q-"max(m,l) ] 


where [|^j implies the least integer K. 

PiRT 0 

COMP MTS ATOR DBSIGR BY A SEQUEROE OP CTO ET^S [?] 


In this the usual way is to assign as many poles as 
possible by a constant feedback matrix E^ (first dyad) and 
then ■fcC' assign the x-emaining poles by second dyad Ep while 
retaining as many as possible poles assigned by feedback E^ . 

)So the desired compensator can be defined as 

E = E^ + Eg ( 01 ) 

where E^ and E^ both are dyadic. 

E^ is determined to place max(m,l) poles at desired 
location as already discussed in part B. Then, the remaining 
problem is to retain as many as possible of tne max(m,l) 
poles, already assigned by the first dyad while attempting to 
shift some (or all) of the remaining [n~max(m,l)] poles to 
the desired locations. 

Oonsider the case of 1 > m. Then using a feedback 


( 02 ) 


A”1 2 


let the resulting Ixm system with 1 poles at the desired 
locations, be expressed as 

/Vhs) (03) 

( 2 ) 

Vife must now' choose ^ such that as many as possible of 
these 1 poles, say h in number, become unobsearvable, that is, 

= £,^(s) ir(2)(e) (04) 

¥ow let us find out what can be the maxiimum value of k. 

(A) 

The choice of f^^ ^ is facilitated by the fact that for 
s equal to any of the roots of A^^^(s) and in particular for 
s equal to any of the roots of (b) ^ the matrix p^^^(s) 
has rank equal -to one. 

So equation (04) leads •’"C a set of k =. m8,x(m,l) equations 
in = max(m,l) unknown parameters of f^^ ' ; 

f (2) ) =0; i = 1, max(m,l) . (05) ? 

Oil [ 

i; 

Taking only one row 'from each for i = 1, i 

max(m,l), eqn. (05) can equally be written as I. 

f^(2) ^ ^ Q (Qgx i 

where ith row of f" corresponds to any non-zero row of 

i 

f 

ft' 

i; 

[' 

I' 

. i, 

f. 
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Now dimension of f is max(m,l) ^ inax(m,l). Eqn. (C 6 ) 
lias a non-^irivial solution only if r lias more rows than 
columns* Tiiis can be acbieved only by assigning a value to 
one (or more) parameters of Thus, we can only restrain 

at the most [max(m,l)-l] poles out of the max(m,l) poles 
assigned hy feedback 

Having determined f^^^^ as above, we have ^ 2 ( 3 ) with 
dimension/ 1 xmin(m,l), with elements having maximum degree 
[(n- 1 )-(max(m,l)- 1 )] = [n-max(m,l) ] . Equation (B4) will now 
become 


X I - 
o o 


^0 


(C7) 


where 6^ has dimension? min(ffi,n~l) *1 , X has dimension 

min(m,n-l +1 ) .m and P = Has dimension mx 1 , 

0 o 

Now if n m-i-l ~1 , then n-1+1 > m. So, the number of 

fjj 

columns of X^ will be greater than or equal^its number of 
rows, and (07) will have a solution given by 




^ -..1 ^ 

= [X J ^ 6 
^0 0 




0 

if 

n = m+1-1 

and 

'P; 

P 

'A 

= [x'^ X 

^ t'K ^ 

^t 

6 ., 

if n < m+l-*1 


But if n> m+1-1 , then eqn. (07) will not have a solution 
unless 6 is a "special case such that it lies in the column 
space of X^ . In general, a dynamic compensator will be needed. 



ip j) lying eqa. (A5) shows that compensator order r 
will he bounded as 


rfeCff, §■??,( ^ ? ll -ljj "-min ( m , 1 ) ^ 

^ max(m,l) J ^ 


< 


[^n-(max(m,l)^ 
min(m,l) ] 


or 


r n— ( m+l-" 1 ) 
^ max { m , 1 ) 


] ^ r 


$ [n~(m+l-»l)] 


for placement of all (n+r) poles of closed loop system. 

Here [Kj implies the least integer K. 

It can also be verified that for placement of (q+r) poles, 
where n $q< mH-l~l ,compensator order r xfill be bounded as 


where [*] denotes least integer but greater than or equal to 
the enclosed quantity. 


PART D 

PULI RANK COMPENSATOR DESIGN BY MINIMAD SEQUENCE 

OE DY/iD.S [10] 

The procedure to be described in this part is essentially 
a generalization of two dyad feedback technique. 

The desired feedback matrix E is constructed from a 


minimal sequence of dyads, that is. 
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F = f ( i) ^ ( i) 

i=t 


(D1) 


where is an mxl vector, is a 1x1 vector, and 

[i = inin(iii,l). The matrix F csn equally be expressed as a 
product of two matrices as 


F=F^Fo, (D2) 

where the columns of F are the vectors f and the rows of 

c c 

Fq are the vectors i = 1»2, ..., p in both cases. The 

matrix F will have maximum rank u, if F and F„ both will 

c o 

have maximum rank. 

In the following discussion, it will be seen that for 

the case m 1, F will have maximum rank by construction. 

c 

It then remains to show that F^ will also have maximum rank. 

The result for the ease of m ,> 1 is obtained by duality. 

Essentially , at each stage in the proposed procedure, 
the vectors are chosen so that one more input to the 

original system is used. It is this action that guarantees 

(i) 

the rank of F^ to be maximum. The vectors f„' ' are then 
c 0 

chosen to assign as many poles as possible at each stage. 

The procedure and its proof is given below. 

let f =[l»0, ... 0]^, say. Then, since it is 
o 

assumed that B has rank m, the pseudo-input vector ^ ^ = Bf ^ ^ ^ 

c 

will influence at least one mode of the system [A,b^^^,C]. let 
,* 1 4n^^n,‘ be the number of modes influenced 
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by b^^^^ The system - has 1 outputs and 1 

input. So now at the most poles can be assigned to 
specific locations at this stage, where q = min(n,,l). 

/ % I I 

Vector fg ^ ^ ' is then determined to place these poles at 
desired locations by solving 

=/\^°^(s) + [adj(sI-A)J b^^^ (1)3) 

where is as defined by eqn. A^°^(s) is the 

open loop characteristic polynomial 0^^^ =fQ^''^0, 

and b^ ^ ^ = Bf ^ ^ \ 
c 

Now let , ..., Oj"^, say, where 3'' 

chosen such that '.as many as possible of the q^ poles assigned 

in stage 1 are made uncontrollable with respect to system 

where A^ ^ ^ = A -■ Bf f and such that 

o o 

/ Q "S 

f^'' ' ^ 0. This ensures that at least one mode of the system 
is controllable, and at least one mode is un- 
controllable. let n^f such that 1 n^ ^ n-1 be the number 
of modes influenced by b^ ^ and let q _2 = min(l,n 2 .'' be the 

number of poles which can be further assigned to specific 

( 2 ) 

values at this stage. The vector f^^ ' is then determined 
to place q 2 poles at desired locations by solving 

= ^(s) + c^^^[adj(sI-A + b^^^c^^^)] b^^^ (D4) 

This procedure is repeated at each subsequent stage until 
at the mth stage the result of second part of Part C of this 
Appendix is applied. 
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It is easy to see that the matrix P = f 

c c 

i = 1,2, n ; generated in this manner jh^has full rank 

equal to m, The fact that is also of full rank has 
been proved by Munro and Hirbod [l0] and hence has been 
omitted here. 

The procedure described here is the basis for the 
algorithm presented in Chapter 4. 
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iPP_.xTOIX B 
MlceiUE MOD El [9] 


A three winding model of synchronous machine is 
considered. The stator of the machine is represented by a 
dependent current source. The power system network is 
represented by a single phase equivalent circuit as shown 
in Eig. 4.2, 


The system non-linear equations are given below : 


2H O ^ 

^0 dt^ “ ® 


dl 

dT 


d 


f^fd - ^d^d ” ^do) 


and 


= e 1 . 


where 


T, 


= X’ (I,i - I i, ) 
d '• d q q d'^ 


and 


(X^ - \)/^d 


(E1) 

(E2) 

(E3) 


The non-state variables Iq expressed 

in terms of 6 and as follows- From Eig. 2.2, the armature 
current phasor 'x can be solved for, by using the network 
equations 


^ q E 

^ R+3TX+X^ ■” R+j(X+Xp 


(a. + ja^)! - (y. + jyp)E 
^ ^ ^ (B4) 
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where 


X^(Z+X^) 

'■■‘■’o'""""'' '' ' " W p 

' R^+(X+Xp"^ 


RX. 


= 


d 


2 E^+(X+X^)^ 


R 

y , = -rr and y^ = -5 ^ P 

1 R^+(X+X|)^ ^ 


R^+(X+Xp 


So from (E 4 ) and Xi^. 2.2 

\ '■ ^^d (2^i + 3a2)(lq~3l(i) " (y^ + 0y2)® ® 

By equating real and imaginary parts 


*06 


(S5) 


iq = ^1^q ^2^d ’ ^2 

and 

id = a^I^ - a^Iq^ + ECy^ cos 6 -• y^ sin 6 ) 


Substituting I from eqns. (E5) in the above equation 

u. 


get 


a. 


B 


■q ' 1 -a.e d l-a.e 


(y^ cos6 + sin6) 


(E 6 ) 


and 


a . I j 


a-gO 


^d ■" I d 1 -a^e 


[a^I^ - E(y^ cos6 + sin6)] 


+ ECyp cos6 - y^ Sind) 


(E 7 ) 


Excitation System Equations 

From the block diagram of excitation system given in 
Eig. 4.3* ^6 state equation can be written as 


we 
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fit 

dt 

. _ -11. ^ _R y 

’’e t 

(E8) 

\d 

” "^ref ' ^‘1 , ^ 

(E9) 


after n^lecting the small time-constant T of static excitation. 

© 

In eqns. (IS) and (E9) 




2A 


't “ ^ ~ ^^'^q ^q^ (E10) 

and u is the output of PSS. 

So the non-linear differential equations can he written as 


u = [s 


fd 


X.'l^ 
' ad 




6 = 


W 


w 


— [T 
2H ^ m 






X 


" 


1 




m + m 

-^R 


(E11 ) 
(E12) 

(E13) 

(B14) 


alongwith 


®fd 


V 


ref 


+ u 


and ig^ and i^ are already defined explicitly, in terms of 
and 6 in ( E6 ) and ( E7 ) . 

Linearized Model 

The linearized model can he given as 
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kd I 

{ } 

! 1 

^12 

0 

Q 

K^^dj 

ii“-’ 

! b 

= 

i 0 

} 

0 

1 

0 j 

( 

i-h6 ; 

+ 1 0 

law i 

I j 

^31 

®32 

0 

0 ; 

ii i 

: i 

/'.w ! 

r 1 

1 

j 0 

i 

i ! I 

ioX^ 1 . 

; ^41 

^"42 

0 

^44 i 

' f 

\t-xA 
; 1 ! 

•i 

! 0 
( 

«Wf.> . i.t,* >’ 

where 






! 




d ""d' 


d 1 


Y ' T” 

•^d do 


^\U 


(B15) 


‘12 


-^d) 
=<:d ’^ao 


a 


14 


^a ’^do 


W 


31 


~ (i' ^ K^I ^ + K.I,^ - K,i,^ q)X' 

2H '■ qo 3 qo 1 do l do d 


"'32 


^4i 


^/T 

^ (^4^0 - ^2^do ^2 "do ^)^d 




T 


(EIC^ ~ XX- ) + ~'3-^ (RI^ + XK^)] 
^to " ' ^to ^ ^ 


a 


R to 
do 


42 


Kn Y 

^ (E qos6 + RK, XK ) 

R to 


V 


+ -02. („E sin6^ + RK„ + XK . ) ] 


V 


to 


0 


a 


44 


R 

1 


' " ’^do 


and 
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Here and are the constants which have been 

obtained after linearizing equations of i^^ and i^ about the 
operating point. 

So, 

Ai = K. fMj + '1-6 

q -I - d 2 

and 

/j,i^ — E^ '^^d E^ /lx6 
where 

2, = . ^2 = Trijr (^1 - 72 

E3 = and = [~(a2£y^ + ( 1~a^ e)y2)siii<5^ 

1 1 “* 8 , ^6 

+ {^^^2 ~ )y^ )cos6^] 

The set of equations (El 5 ) can be written in the following 
standard form 

m 

X = Ax + bv 
"t/ 

where x = A 6, /uw, hx^] is the 4th order state vector and 

V =i [ A-u] is the scalar input which is the output of PSS too. 

The matrices A and b are constant matrices at 
operating point, and their numerical values are given as 


B 


'"■'-0.16312 

-0.3773 

0.0 

-0.34722" 

0.0 

j 

0.0 

1.0 

0.0 

-9.0999 

-25.7124 

0.0 

0.0 

1983.2893 - 

■2436.8389 

0.0 

-100.0 


0. 34722"] 


b 


0.0 

0.0 



System Data 


Generator Prameters ; 


= 1.14, X^ = 0.66, X^ = 0.24 

^do = H = 1 .5, Wq = 314.0 


Uetwork Paraviieters s 

R = 0.0, X = 0.7 


Operating Data % 

S = 1 *0 Zo° 

P = Received power = 1.0 
Q = Received reactive power =0.0 

= Voltage regulator time constant=0.0l sec. 



